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ABSTRACT 

We report on a study of the parameters characterizing the mass and velocity 
distributions of two samples of relaxed elliptical-like-objects (ELOs) identified, 
at z = 0, in a set of self-consistent hydrodynamical simulations operating in 
the context of a concordance cosmological model. Star formation (SF) has been 
implemented in the simulations in the framework of the turbulent sequential 
scenario through a phenomenological parameterization that takes into account 
stellar physics processes implicitly through the values of a threshold gas density 
and an efficiency parameter. Each ELO sample is characterized by the values 
these parameters take. We have found that the (logarithms of the) ELO stel- 
lar masses, projected half- stellar mass radii, and stellar central l.o.s. velocity 
dispersions define dynamical Fundamental Planes (FPs). Zero-points depend on 
the particular values that the SF parameters take, while slopes do not change. 
The ELO samples have been found to show systematic trends with the mass 
scale in both, the relative content and the relative distributions of the baryonic 
and the dark mass ELO components. The physical origin of these trends lies in 
the systematic decrease, with increasing ELO mass, of the relative dissipation 
experienced by the baryonic mass component along ELO mass assembly, result- 
ing into a tilt of the dynamical FP relative to the virial plane. ELOs also show 
kinematical segregation, but it does not appreciably change with the mass scale. 
We have found that the dynamical FPs shown by the two ELO samples are con- 
sistent with that shown by the SDSS elliptical sample in the same variables, with 
no further need for any relevant contribution from stellar population effects to 
explain the observed tilt. These effects could, however, have contributed to the 
scatter of the observed FP, as the dynamical FPs have been found to be thinner 
than the observed one. The results we report on hint, for the first time, to a 
possible way to understand the tilt of the observed FP in a cosmological context. 

Subject headings: dark matter: galaxies: elliptical and lenticular, cD - galaxies: 
evolution - galaxies: formation - galaxies: fundamental parameters - hydrody- 
namics - methods: n-body simulations 
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1. INTRODUCTION 

The 3-parameter space of the observed effective radius, i?g lght , the mean surface bright- 
ness within that radius, < I hght > e , and the central line-of-sight (los) velocity dispersion, 
cios.O) °f early-type galaxies is not homogeneously populated. These galaxies define a plane 
on this space, known as the Fundamental Plane (FP, Djorgovski & Davis 1987; Dressier et 
al. 1987a; Faber et al. 1987; Kormendy & Djorgovski 1989) defined by: 

log 10 R^ ht = a log 10 a los , + b log 10 < J light > e +c. (1) 

A new standard of reference for nearby elliptical galaxies is provided by the Sloan digital 
sky survey (SDSS, see York et al. 2000) sample of early-type galaxies (see Bernardi et al. 
2003a, 2003b, 2003c for the Early Data Release), containing to date 9000 morphologically 
selected ellipticals from different environments, a number larger than the number of ellipticals 
in all the previously analyzed samples. The values of the FP coefficients from this sample 
are a ~ 1.5, similar in the four SDSS bands, b ~ —0.77, and c ~ —8.7 (see their exact values 
in Bernardi et al. 2003c, Table 2) with a small scatter. These SDSS results confirm previous 
ones, either in the optical (Lucey, Bower & Ellis 1991; de Carvalho & Djorgovski 1992; 
Bender, Burstein & Faber 1992; Jorgensen et al. 1993; Prugniel & Simien 1996; Jorgensen 
et al. 1996) or in the near-IR wavelenghts (Recillas-Cruz et al. 1990, 1991; Pahre, Djorgovski 
& de Carvalho 1995; Mobasher et al. 1999), even if the published values of a show larger 
values in the iCband than at shorter wavelengths (see, for example, Pahre, de Carvalho & 
Djorgovski 1998). 

The existence of the FP and its small scatter has the important implication that it 
provides us with a strong constraint when studying elliptical galaxy formation and evolution 
(Bender, Burstein & Faber 1993; Guzman, Lucey & Bower 1993; Renzini & Ciotti 1993). 
The physical origin of the FP is not yet clear, but it must be a consequence of the physical 
processes responsible for galaxy assembly. These processes built up early type galaxies as 
dynamically hot systems whose configuration in phase space are close to equilibrium. Taking 
an elliptical galaxy as a system in equilibrium, a mass scale Mfp can be estimated from the 
observable i^ ight , < J light > e and cios,o parameters through the expression: 

M FP = 3c M al sfi R^ u /G, (2) 
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where G is the gravitational constant and cm is a mass structure coefficient. The total 
luminosity of the galaxy is given by the expression L = 2vr < I li % ht > e (i? u § ht ) 2 . As Faber et 
al. (1987) first pointed out, the virial theorem M vir = cp^ *) 2 ?-* * /G, with a^ ot the average 
3-dimensional velocity dispersion of the whole elliptical, including both dark and baryonic 
matter, r* ot the dynamical half-radius or radius enclosing half the total mass of the system, 
and Cf a form factor of order unity, would imply 

plight _ 3c M r -k°fos,0 < -f'' ght > e 1 , s 

Ke ~ 2nGM viv [6) 

where c^ r is defined by Eq. (2) when M F p = M vir , that is 



vir _ GM vn 



or 

C M = c fc v c t (5) 

with c v = (^3 0t ) 2 /3of os ,o an d c r — ^e 0t /-^e ght - Assuming that both the quantities and the 
dynamical mass-to-light ratios, M vir /L, are independent of mass, the scaling relation i? hght oc 
°foso < J llght >~ l would then hold. However, this predicted scaling law is inconsistent with 
those found observationally (a ^ 2,6 ^ — 1), i.e., the FP is tilted relative to the virial 
relation, implying that at least some of the assumptions made to derive it is incorrect. 

Different authors interpret the tilt of the FP relative to the virial relation as caused 
by different misassumptions that we comment briefly (note that we can write M v ; r /L = 
M star /L x M vir /M star , where M star is the stellar mass of the elliptical galaxy): i) A first 
possibility is that the tilt is due to systematic changes of stellar age and metallicity with 
galaxy mass, or, even, to changes of the slope of the stellar initial mass function (hereafter, 
IMF) with galaxy mass, resulting in systematic changes in the ste//ar-mass-to-light ratios, 
M star /L, with mass or luminosity (Zepf & Silk 1996; Pahre at al. 1998; Mobasher et al. 
1999). But these effects could explain at most only ~ one third of the (3 ^ value in the 
-B-band (Tinsley 1978; Dressier et al. 1987; Prugniel & Simien 1996; see also Renzini & 
Ciotti 1993; Trujillo, Burkert & Bell 2004). Futhermore, early-type galaxies in the SDSS 
have been found to have roughly constant stellar-mass-to-light ratios (Kauffmann et al. 
2003a, 2003b). Anyhow, the presence of a tilt in the i^-band FP, where population effects 
are no important, indicates that it is very difficult that the tilt is caused by stellar physics 
processes alone, as Bender et al. (1992), Renzini & Ciotti (1993), Guzman et al. (1993), 
Pahre et al. (1998), among other authors, have suggested, ii) A second possibility is that 
M vil /L changes systematically with the mass scale because the total dark-to- visible mass 
ratio, M vir /M star changes (see, for example, Renzini & Ciotti 1993; Pahre et al. 1998; 
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Ciotti, Lanzoni & Renzini 1996; Padmanabhan et al. 2004). Otherwise, a dependence of 
on the mass scale could be caused by systematic differences in iii), the dark versus 
bright matter spatial distribution, which could be measured through systematic variations 
of the c r coefficients with mass, iv), the kinematical segregation, the rotational support 
and/or velocity dispersion anisotropy in the stellar component (dynamical non- homology), 
measurable through the cy coefficients, and, v), systematic geometrical effects, measurable 
through the c F , c r or c v coefficients. Taking into account these effects in the FP tilt demands 
modellizing the galaxy mass and velocity three-dimensional distributions and comparing the 
outputs with high quality data. Bender et al. (1992) considered effects iii) and iv); Ciotti 
et al. (1996) explore ii) - iv) and conclude that an systematic increase in the dark matter 
content with mass, or differences in its distribution, as well as a dependence of the Sersic 
(1968) shape parameter for the luminosity profiles with mass, may by themselves formally 
produce the tilt; Padmanabhan et al. (2004) find evidence of effect ii) in SDSS data. Other 
authors have also shown that allowing for broken homology, either dynamical (Busarello et 
al. 1997), in the luminosity profiles (Trujillo et al. 2004), or both (Prugniel & Simien 1997; 
Graham & Colless 1997; Pahre et al. 1998), brings the observed FP closer to the virial plane. 

One important source of ambiguity in observational data analysis comes from the impos- 
sibility to get accurate measurements of the elliptical three-dimensional mass distributions 
(either dark, stellar or gaseous) and velocity distributions. Analytical models give very inter- 
esting insights into these distributions as well as the physical processes causing them, but are 
somewhat limited by symmetry considerations and other necessary simplifying hypotheses. 
Self-consistent gravo-hydrodynamical simulations are a very convenient tool to work out this 
problem, as they directly provide with complete 6-dimensional phase-space information on 
each constituent particle sampling a given galaxy-like object formed in the simulation, that 
is, they give directly the mass and velocity distributions of dark matter, gas and stars of 
each objet. This phase space information allows us to test whether or not the c^ r (that is, 
the cf, and c r ) coefficients, as well as the M vir /M star ratios, do or do not systematically de- 
pend on the mass scale. This is the issue addressed in this Letter, where we analyze whether 
the dependence is such that the tilt and the scatter of the observed FP can be explained 
in terms of the regularities in the structural and dynamical properties of ELOs formed in 
self-consistent hydrodynamical simulations. 

2. The Dynamical Fundamental Plane of Simulated Objects 

In a self-consistent numerical approach initial conditions are set at high z as a Monte- 
carlo realization of the field of primordial fluctuations in a given cosmological model; then 
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the evolution of these fluctuations is numerically followed up to z = by means of a com- 
puting code that solves the N-body plus hydrodynamical evolution equations. We have run 
ten simulations in the framework of a flat ACDM cosmological model, with Q\ = 0.65, 
^baryon = 0.06, o" 8 = 1.18 and h = 0.65. The code used in our simulations is DEVA (Serna, 
Dommguez-Tenreiro, & Saiz 2003). In this code, particular attention has been paid to the 
implementation of conservation laws (energy; entropy, taking into account the V/i terms; an- 
gular momentum). Star formation (SF) has been implemented in the code in the framework 
of the turbulent sequential scenario (Elmegreen 2002) through a phenomenological parame- 
terization, that transforms cold locally-collapsing gas, denser than a threshold density, p t hrcs, 
into stars with a timescale given by the empirical Kennicutt-Schmidt law (Kennicutt 1998), 
with an average star formation efficiency at the scales resolved by the code c*; possible feed- 
back effects are implicitely taken into account through the values of the SF parameters. In 
any run, 64 3 dark matter and 64 3 baryon particles, with a mass of 1.29 x 10 8 and 2.67x 10 7 M Q , 
respectively, have been used to homogeneously sample the density field in a periodic box of 
10 Mpc side. We refer the reader to Serna et al. 2003 and to Saiz, Dommguez-Tenreiro, 
& Serna 2004 for further details on the simulation technique and SF implementation in the 
code. Five out of the ten simulations (the SF-A type simulations) share the SF parameters 
(pthrcs = 6 x 10~ 25 gr cm~ 3 , c* = 0.3) and differ in the seed used to build up the initial 
conditions. To test the role of SF parameterization, the same initial conditions have been 
run with different SF parameters (pthres = 1-8 x 10 -24 gr cm -3 , c* = 0.1) making SF more 
difficult, contributing another set of five simulations (hereafter, the SF-B type simulations). 
Galaxy-like objects of different morphologies appear in the simulations. ELOs have been 
identified as those objects having a prominent dynamically relaxed stellar spheroidal compo- 
nent, with no disks and very low cold gas content. This stellar component has typical sizes 
of no more than ~ 10-40 kpc and it is embedded in a halo of dark matter typically ten 
times larger in size. ELOs have also an extended corona of hot diffuse gas. It turns out that 
26 (17) out of the more massive objects formed in SF-A (SF-B) type simulations fulfil this 
condition, giving the hereafter termed SF-A and SF-B ELO samples. Note that due to their 
respective SF implementations, galaxy-like objects formed in SF-A type simulations tend to 
be of earlier type than their counterparts formed in SF-B type simulations; this is why the 
ELO number in SF-A sample in higher than in SF-B sample. Moreover, gas has had more 
time to lose energy along SF-B type ELO assembly than in their SF-A type counterparts, 
and, consequently, the former have smaller sizes than the latter. 

Following the discussion in §1, two mass scales have been measured on ELOs: the virial 
mass M vir , or total mass of the ELO at its halo scale (we adopt the fitting formula of Bryan 
& Norman, 1998, for the spherical overdensity at virialization), and the stellar mass at the 
ELO scale, M^* ar . Concerning length scales, the relevant ones are: i), at the halo scale: the 



- 6- 



virial radii, r vir and the half-total mass radii, r*°^, enclosing M vir /2; and ii), at the ELO scale, 
the stellar half-mass radii, rf^ Q , defined as those radii enclosing half the M^* ar mass, and 
the projected stellar half-mass radii, Rf^ T Q , measured onto the projected mass distribution. 
The total (including baryons and dark matter) mean square velocity within r v ; r , cr^, as 
well as the mean square stellar and central stellar 2 l.o.s velocity dispersions, (rf™ and of*^, 
respectively, have also been measured on ELOs. Note that the scales entering the virial 
relation are and 0-3°^ and that they are not observationally available. Assuming that 
the projected stellar mass distribution, S star (i?), can be taken as a measure of the surface 
brightness profile, then < S star >c= c < J light > e , with c a constant, and Rf^ Q ~ i?J, ight 
and we can look for a fundamental plane (hereafter, the dynamical FP) in the 3-space of 
the structural and dynamical parameters Rf^ r Q , < S star > c and crf^ , directly provided by 
the hydrodynamical simulations. To make this analysis as clear as possible, we transform 
to a K-like orthogonal coordinate system, the dynamical nf system, i— 1,2,3, similar to that 
introduced by Bender, Burstein & Faber (1992), but using Rf^ r Q instead of i?j! ght and < 
Sstar > c instead of < J llght > e , and, consequently, free of age, metallicity or IMF effects. The 
k d and k coordinates are related by the expressions: Ki ~ k^, k 2 ~ k% — a/6/3 logM^* ar /L 
and k 3 ~ + a/3/3 log M^ r /L. We discuss the tilt and the scatter of the dynamical FP 
separately. We first address the tilt issue. We use at this stage for the Rf^o an< ^ a i S os,o 
variables the averages over three orthogonal l.o.s. projections, to minimize the scatter in the 
plots caused by projection effects. 

Figure 1 plots the versus «° (top) and k% versus /tf (bottom) diagrams for ELOs 
in the SF-A (filled symbols) and SF-B (open symbols) samples. We also drew the 2a con- 
centration ellipses in the respective variables, as well as its major and minor axes, for the 
SDSS early-type galaxy sample in the z band as analyzed by Bernardi et al. 2003b, 2003c 3 . 
We recall that the ellipse major axis corresponds to the orthogonal mean square regression 
line for the two variables in the Figure (for further details see Saiz et al. 2004). The most 
outstanding feature of Figure 1 (upper panel) is the good scaling behaviour of versus /tf , 
with a very low scatter (see the slopes Mi in Table 1; note that the slopes for the SF-A 
and SF-B samples are consistent within their errors while the zero-points depend on the SF 
parameterization through the ELO sizes). The values of the slopes in Table 1 mean that 
systematic variations of the structural and dynamical properties of ELOs with the mass 
scale cause, by themselves, a tilt of the dynamical FP relative to the virial relation. Another 
interesting feature of Figure 1 is that it shows that most of the values of the nf coefficients 



2 Rccall that the empirical l.o.s velocity dispersion, (Ti os ,o, is measured through stellar spectra 

3 The constant stellar-mass-to-light ratios allow us to write the covariance matrix using the E = logM^ ar 
variable instead of absolute magnitude or the logarithm of the luminosity L 
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are whithin the 2a concentration ellipses in both plots for ELOs formed in SF-A type simula- 
tions, with a slightly worse agreement for ELOs in the SF-B sample. This means that ELOs 
have counterparts in the real world (Saiz et al. 2004). It is also worth mentioning that these 
results are stable against slight changes in the values of the Q\, ^baryon or h parameters; for 
example, we have tested that using their preferred WMAP values shows results negligibly 
different to those plotted in Figure 1. 

To deepen into the causes of the tilt, we have calculated the slope f3 v - ir of the M vir /M^* ar oc 
(M b i * ar ) /3vir scaling relation for ELOs in the SF-A and SF-B samples (Table 1). We got 
/3 v i r > 0, indicating that the mass fraction of stars bound to the ELOs (or, more gener- 
ically, cold baryons) relative to total mass within r vir decreases with the mass scale (as 
suggested by Renzini & Ciotti 1993 and Pahre et al. 1998). We have also found that 
the c^ r coefficients show a mass dependence that can be parametrized as a scaling relation 
c^ r oc (M^)^ (broken homology, see Table 1). As discussed in §1, different possible sources 
for homology breaking exist. To dilucidate which of them are relevant, we first note that 
the cy and c r coefficients can be written as: cy = cvd x cypc, with cvd = (°3°h/ ' a fbo) 2 an d 
cvpc = « b a o) 2 /3«^) 2 and c r = c rD x c rP with c rD = r™/r% and c rP = r£/H*&. This 
gives cj^ r = cpCvDCvpc c rDC r p ; where the cvd and c r p factors measure the kinematical and 
spatial segregations between dark and stellar matter, respectively, while cypc an d c r p mea- 
sure projection and other geometrical effects in the stellar mass and velocity distributions. 
Writting Ci oc (Mjf G ar ) ft , with i = F, VD, VPC, rD, rP, we get (3 M = /3f+/3vd+/3vpc+&d+&p 
when the $ slopes are calculated through direct fits. These $ slopes are given in Table 1, 
as well as their 95 % confidence intervals both for the SF-A and SF-B samples. We see that, 
irrespective of the SF parameterization, the main contribution to the homology breaking 
comes from the c r p coefficients (see Guzman et al. 1993), while (3p and (3yd have values 
consistent with cp and cvd being independent of the ELO mass scale, c r p and cypc show a 
very mild mass dependence in the SF-A sample and none in the SF-B sample. 

We now turn to consider the scatter of the dynamical FP for the ELO samples and 
compare it with the scatter of the FP for the SDSS elliptical sample, calculated as the 
square root of the smallest eigenvalue of the 3x3 covariance matrix in the E (or logL), V 
and R variables (Saglia et al. 2001). As Figure 1 (top) suggests, when projection effects 
are circunvented by taking averages over different directions, the resulting three dimensional 
orthogonal scatter for ELOs is smaller than for SDSS ellipticals (opvr = 0.0164 and <tevr = 
0.0167 for the SF-A and SF-B samples, respectively, to be compared with olvr = 0.0489 
for the SDSS in the logL, V and R variables). To estimate the contribution of projection 
effects to the observed scatter, we have calculated the orthogonal scatter for ELOs when no 
averages over projection directions for the Rf^l an d ^loTo variables are made. The scatter 
(o"evr = 0.0238 and opvr = 0.0214 for the SF-A and SF-B samples) increases, but it is still 
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lower than observed. This indicates that a contribution from stellar population effects is 
needed to explain the scatter of the observed FP, as suggested by different authors (see, for 
example, Pahre et al. 1998; Trujillo et al. 2004). 

To sum up, the ELO samples have been found to show systematic trends with the mass 
scale in both, the relative content and the relative distributions of the baryonic and the 
dark mass ELO components. These trends do not significantly depend on the star formation 
parameterization and they are due to a systematic decrease, with increasing ELO mass, of 
the relative amount of dissipation experienced by the baryonic mass component along ELO 
formation, a possibility that Bender et al. (1992), Guzman et al. (1993), and Ciotti et al. 
(1996) had suggested. These trends cause a tilt of the virial plane in such a way that there 
is no further need of any relevant contribution from stellar population effects to explain the 
observed tilt. The scatter of the observed FP, however, probably requires a contribution 
from such stellar effects. 
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Fig. I. — Edge-on projection (top panel) and nearly-face-on projection (bottom panel) of the 
dynamical FP of ELOs in the k d variables (filled circles: SF-A sample; open circles: SF-B 
sample). We also draw the respective concentration ellipses (with their major and minor 
axes) for the SDSS early-type galaxy sample from Bernardi et al. (2003b) in the z-band. 
See text for more details 
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Table 1: 





SF-A 




SF-B 




Mi 


0.256 


± 0.035 


0.281 


± 0.048 




0.221 


± 0.083 


0.237 


± 0.158 


/3m 


-0.204 


± 0.116 


-0.247 


± 0.189 




0.025 


± 0.048 


0.022 


± 0.081 




0.021 


± 0.041 


0.076 


± 0.075 


Pvpc 


-0.044 


± 0.029 


-0.044 


± 0.093 




-0.225 


± 0.127 


-0.316 


± 0.199 


PrP 


0.019 


± 0.009 


0.016 


± 0.017 



Column 2: the slopes of the k° = M 1 k^+M relation (direct fits); the slopes of the M vir /M^* ar 
and q oc (M^* ar ) ft scaling relations for the the SF-A sample, calculated in log — log plots 
through direct fits. Column 3: their respective 95% confidence intervals. Columns 4 and 5: 
same as columns 2 and 3 for the SF-B sample. 



